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ABSTRACT 

Weak gravitational lensing induces distortions on the images of background galaxies, 
and thus provides a direct measure of mass fluctuations in the universe. The distortion 
signature from large-scale structure has recently been detected by several groups for 
the first time, opening promising prospects for the near future. Since the distortions 
induced by lensing on the images of background galaxies are only of a few percent, 
a reliable measurement demands very accurate galaxy shape estimation and a care- 
ful treatment of systematic effects. Here, we present a study of a shear measurement 
method using detailed simulations of artificial images. The images are produced us- 
ing realisations of a galaxy ensemble drawn from the HST Groth strip. We consider 
realistic observational effects including atmospheric seeing, PSF anisotropy and pix- 
elisation, incorporated in a manner to reproduce actual observations with the William 
Herschel Telescope. By applying an artificial shear to the simulated images, we test the 
shear measurement method proposed by Kaiser, Squires & Broadhurst (1995, KSB). 
Overall, we find the KSB method to be reliable with the following provisos. First, 
although the recovered shear is linearly related to the input shear, we find a coeffi- 
cient of proportionality of about 0.8. In addition, we find a residual anti-correlation 
between the PSF ellipticity and the corrected ellipticities of faint galaxies. To guide 
future weak lensing surveys, we study how seeing size, exposure time and pixelisation 
affect the sensitivity to shear. We find that worsened seeing linearly increases the noise 
in the shear estimate, while the sensitivity depends only weakly on exposure time. The 
noise is dramatically increased if the pixel scale is larger than that of the seeing. In 
addition, we study the impact of overlapping isophotes of neighboring galaxies, and 
find that this effect can produce spurious lensing signals on small scales. We discuss 
the prospects of using the KSB method for future, more sensitive, surveys. Numerical 
simulations of this kind are a required component of present and future analyses of 
weak lensing surveys. 

Key words: cosmology: observations, gravitational lensing, large-scale structure of 
Universe, methods: data analysis, techniques: image processing 



1 INTRODUCTION 

Weak lensing provides a unique method to directly measure 
the mass fluctuations on large scales in the universe (see 
Melher 1999; Kaiser 1999; Bartelmann & Schneider 1999 for 
recent reviews). This method relies on the measurement of 
small, coherent distortions produced by lensing upon the 
shapes of background galaxies. This effect is now routinely 
used to map the mass of clusters of galaxies (see reviews 
by Fort & MeUier 1994, Schneider 1996). Recently, the tech- 
nique was extended to the field by several groups who re- 
ported the statistical detection of weak lensing by large- 
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scale structure (Wittman et al. 2000; van Waerbeke et al. 
2000; Bacon, Refregier & Ellis 2000 (ERE); Kaiser, Wil- 
son & Luppino 2000). More precise measurements of this 
"cosmic shear" from upcoming observations will provide in- 
valuable cosmological information (eg. Kaiser 1992; Jain & 
Seljak 1997; Kamionkowski et al. 1997; Kaiser 1998; Hu & 
Tegmark 1998; van Waerbeke et al. 1998). 

Because the distortions induced by lensing are only of 
the order of 1%, these measurements are very challenging. 
In particular, they require tight control of systematic effects 
and a precise method for the measurement of the shear. One 
of the potential weaknesses of the cosmic shear programme 
is the step leading from the measurement of the shapes of 
galaxies to the estimation of the lensing shear, in the pres- 
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ence of an anisotropic PSF. The first metfiod proposed to 
treat this problem was that by Bonnet & Mellier (f995). A 
more general, and now widely used, method was proposed 
by Kaiser, Squires & Broadhurst (f995, KSB) and further 
developed by Luppino & Kaiser (1997) and Hoekstra et al. 
(1998). Variations and alternatives to the KSB method have 
since been presented by Kaiser (1999b), Rhodes et al. (1999) 
and Kuijken (1999). 

In this paper, we address the ellipticity-to-shear prob- 
lem using numerical simulations of artificial images. The 
numerical simulations which have been used for this pur- 
pose in the past have been derived either from HST images, 
degraded to match ground-based observations (eg. KSB, 
Wittman et al. 2000), or by using ab-initio artificial galaxy 
catalogs (eg. Kaiser 1999b). The former approach provides 
accurate shape statistics for the simulated galaxies, but can 
only produce a small simulated area. The latter approach al- 
lows the simulation of arbitrarily large areas, but is not nec- 
essarily as realistic. Because we are aiming at the demanding 
cosmic shear regime, we thus use a hybrid method in which 
large realisations of artificial galaxy images are drawn to 
reproduce the statistics of existing HST surveys. 

Because it is widely used and more documented, we 
focus on the KSB method, and test its reliability in realis- 
tic observational conditions. For definitiveness, we consider 
the weak-lensing survey of BRE, who used ground based 
observations with the William Herschel Telescope (WHT). 
We produce artificial galaxy catalogs generated from ran- 
dom realisations based on the HST Groth Strip (Groth et 
al. 1994, Rhodes 1999). By applying artificial shears to the 
simulated images, we test both the systematic and statisti- 
cal uncertainties of the method. We also investigate how the 
shear signal degrades as a function of seeing, exposure time 
and pixel size. This is of considerable practical interest for 
the design of future weak lensing surveys. In addition, we 
examine the impact of overlapping isophotes on the shear 
signal, an effect which can potentially limit weak lensing 
measurements on small scales. An independent study of the 
shear measurement method using numerical simulations is 
presented in Erben et al. (2000) 

The paper is organised as follows. In we describe 
how we generate the artificial galaxy and star catalogues. In 
§^ we show how these are used to produced realistic images, 
which are compared to observed WHT images. This is fol- 
lowed in §^ be a brief description of our implementation of 
the KSB method. In f|^, we present our results; the accuracy 
of recovery of shear with the KSB method is discussed, and 
the occurrence of an anti-correlation of shear with star ellip- 
ticity at low level is noted. We demonstrate the degradation 
of signal-to-noise with increasing seeing, exposure time and 
pixel size, and also discuss the level at which overlapping 
isophotes will enhance the cosmological shear signal. The 
results are discussed and summarized in Spl 



2 SIMULATED OBJECT CATALOGUE 

The first step in these simulations is to construct an object 
catalogue. To do so, we used the image statistics from the 
Groth Strip, a deep survey taken with the Hubble Space 
Telescope (Groth et al. 1994, Rhodes et al 1999). This HST 
survey is sampled at 0.1 arcsec and thus effectively gives 



us the unsmeared (i.e. before convolution with ground-level 
seeing) ellipticities and diameters of an ensemble of galaxies 
suitable for simulations of ground-based observations. The 
survey consists of a set of 28 contiguous pointings in V and /, 
with an area of approx. 108 arcmin^ in a 3'. 5 x 44'. region. 
The magnitude limit is 7 ~ 26 (WFPC2 I band, F814W), 
and the strip includes about 10,000 galaxies. 

We use a SExtractor (Berlin & Arnoults 1996) cat- 
alogue derived from the entire strip by Ebbels (1998). 
It contains, for each object, a magnitude determined by 
aperture photometry, and a diameter and ellipticity de- 
rived from second-order moments using a top-hat weight 
function. Armed with this catalogue, we model the multi- 
dimensional probability distribution of galaxy properties (el- 
lipticity - magnitude - diameter) sampled by this catalogue. 
We find that the differential galaxy counts as a function 
of /-magnitude are well described by ^ ~ iq('-'0+^iI+'^2i'^) 
galaxies deg~^ mag~^, with ao — —19.0, ai = 1.64 and 
a2 — —0.027; the radius distribution is modelled as r = 
0.095 X io^bi+b2i+b3) arcsecs with 62 = -0.14 and with bi a 
gaussian distribution with mean 3.75 and rms 0.098. Ellip- 
ticity components ei and 62 are described by gaussian prob- 
ability distributions of rms 0.34, and the position angle is 
randomly chosen. We draw from this modeled distribution 
a catalogue of galaxies statistically identical to the Groth 
strip distribution by Monte Carlo selection. 

For definiteness, we aim to reproduce the conditions 
of our weak lensing survey derived from observations with 
WHT (see BRE). This survey consisted of 14 independent 
blank fields observed the WHT prime focus CCD Camera 
(field of view 8' x 16', pixel size 0.237", EEV CCD) in the 
R band. A relevant issue is the number of stars obtained: by 
tuning the galactic latitude (30° < b < 70°), we required the 
fields to contain ~ 200 stars with 7? < 22 in order to map 
carefully the PSF and the camera distortion across the field 
of view. The integration time on these fields was 1 hour, 
affording a magnitude limit of 7?=25.2 (all R magnitudes 
quoted as Vega magnitudes). 

Given that the Groth strip is in / while our data is in 
R, we allowed a slight increase (multiplication by 1.2) to 
radius with magnitude to better model the WHT images; 
this factor has been included in the described model above. 
The number density-magnitude dependence was found to fit 
very well without alteration. 

We spatially distribute the galaxies with a uniform 
probability across the field of view. Since the Groth strip 
does not contain enough stars to create a good model, 
star counts with 7?-magnitude are modeled as a power law 
(Jl oc R°-^) which is found to be a good fit to the WHT 
data. 

We assign a morphological class (elliptical, spiral or ir- 
regular) for future image realisation to each galaxy using the 
results of Abraham et al (1996). Specifically, we use their 
measured fraction for each class as a function of magnitude. 
After application of all of the above procedure, we obtain 
our unlensed object catalogues. 

To produce the lensed object catalogue, we sheared the 
galaxies in the catalogue by calculating the change in the 
object ellipticity due to lensing. This was done using the 
relation (Rhodes et al 1999): 

e- = ei 4- 2(5,j - eiCj)7j (1) 
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Since we are primarily interested only in the mean shear 
measured on the field, we chose the imposed shear to be 
uniform over a given field. 

Stellar ellipticities (simulating tracking errors, atmo- 
spheric effects, etc) are similarly chosen as uniform over 
a given field. The ellipticity for each field is taken from a 
Gaussian probability distribution with a standard deviation 
of o"e = 0.08. This is conservatively chosen to be slightly 
worse than the rms stellar ellipticity of the stars in our WHT 
survey, for which ~ 0.07. 



3 IMAGE REALISATION 

We create the artificial images using the IRAF Artdata 
package. This takes the star and galaxy catalogues, and 
plots the objects with specified positions, ellipticity, mag- 
nitude, diameter and morphology. Only exponential discs 
and de Vaucouleurs profiles are supported. We model ellip- 
ticals and irregulars as de Vaucouleurs profiles, and spirals 
as exponential discs. 

Each pixel in these simulations is subdivided into a 
10 X 10 grid of subpixels. The appropriate subpixel fiux for 
a star or galaxy is calculated from the analytical intensity 
profile, and the PSF convolution is similarly performed at 
the subpixel level. 

We use the package to recreate several WHT-specific 
details: the magnitude zero point is chosen to match the 
telescope throughput, the stars and galaxies are convolved 
with the chosen elliptical PSF (seeing chosen to be 0.8" 
unless otherwise specified, ellipticity dispersion 0.08), the 
image is appropriately pixellised (0.237" per pixel), Pois- 
son photon noise for objects and sky background are in- 
cluded, and gaussian CCD read noise (3.9 electrons) are 
added. The appropriate gain (1.45 electrons / ADU) is in- 
cluded, and an appropriate sky background (10.7ADU per 
sec) is imposed. The PSF profile chosen is the Moffat pro- 
file, I{r) = {1 + {2^"^ - l)(r/rscaio)^)~'^, where /3 = 2.5 and 
rscaie is the seeing radius. The generalised radius r is the 
distance from the centroid, transformed so that the profile 
is elliptical. This profile has wings which fall off more slowly 
than for a gaussian profile, and provides a good description 
of our seeing-dominated PSF. 

An example 4' x 4' portion of a 16 x 8' simulated field 
is shown in figure ^ It can be compared to an observed 
WHT image shown in figure |^ A global impression can be 
obtained from the full 16 x 8' fields plotted in BRE. The 
simulated image lacks saturated stars [R ^ 18); these are, by 
construction absent from the simulated catalogue (see size- 
magnitude figure in BRE) , since such stars cannot be used in 
our weak lensing analysis, and would be immediately excised 
if present. The galaxy images appear to be very similar in 
the simulations as compared with the data. 

To compare these images quantitatively, we derived a 
measured object catalog using the Imcat routines (see §| be- 
low for details). They provide the position, magnitude, half- 
light radius, ellipticity and polarisability tensors for each 
object detected on the image. 

The resulting distribution on the radius-magnitude 
plane is shown in section 7 of BRE for both a simulated 
field and the observed WHT field; the running mean and 
standard deviation of Tg with magnitude is shown on figure 
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Figure 1. Detail of a real data image (WHT3). The area dis- 
played is 4'x 4', while the full image is 8' x 16'. 
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Figure 2. Detail (4'x 4') of a simulated image. 



^ The distributions are similar, with mean radius agreeing 
to within 0.2 pixels throughout the magnitude range. The 
rms scatter in radius is somewhat larger at faint magnitudes 
for the real data, due to the response of the measurement 
software to the simple smooth profiles used in the simula- 
tions. 

Figure ^ compares differential number counts with R- 
magnitude for the simulations and real data. The counts 
derived from the simulated image and from our real _R-band 
fields agree very well; the simulated image counts also agree 
well with the simulated catalogue counts, with an unsur- 
prising turn-over near the expected magnitude completeness 
limit for the 1 hour exposure time {R ^ 25). The simula- 
tions' close impersonation of real data number counts is of 
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Figure 3. Galaxy size statistics with R magnitude. The upper 
curves show the mean filter radius rg in unit magnitude bins; 
solid line shows real data, dotted line shows simulated data. The 
lower curves show the standard deviation on the radius in unit 
magnitude bins. 
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Figure 4. Differential galaxy number counts with ij-magnitude 
(except for real Groth data in /). Real WHT data (WHT3) is 
shown as solid line; Groth strip number counts are shown as dash- 
dotted. The counts for the simulated catalogue are shown dashed, 
and those recovered by KSB from the simulated images is shown 
dotted. 



importance for obtaining realistic results for, eg, overlapping 
isophotes in section ^. 

Figure ^ shows the ellipticity (ei) distributions /(ei) 
for the initial simulated catalogue, and the smeared simu- 
lated and real objects. The distributions were normalised so 
that J" dei /(ei) = 1. As expected, smearing reduces the 
ellipticity dispersion. The simulated and real smeared dis- 
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Figure 5. Comparison of shape measures. The top panel shows 
the normalised ellipticity e\ distribution for initial unsmeared 
simulated catalogue (solid) and smeared simulated (dotted) and 
real (dashed) objects. The bottom panel shows the distribution 
of the shear estimators for corrected simulated (dotted) and real 
(dashed) objects. 



tributions are remarkably similar. See section ^ for further 
discussion on this point. 



4 SHEAR MEASUREMENT METHOD 

Our shear measurement method is a version of the KSB 
method and was described in detail in BRE. Here we sum- 
marize it and then describe the relevant details of our specific 
implement at ion . 



4.1 Overview of the KSB method 

The KSB method derives the shear from the ellipticity of 
galaxies, after correcting for the smearing by the PSF. We 
use the implementation of the KSB method achieved by the 
imcat software kindly provided to us by Nick Kaiser. 

The routine hf indpeaks first finds objects in each field 
by convolving the image with smoothing kernels of different 
sizes. The object radius Vg is defined by the size of the kernel 
which maximizes the signal-to-noise v of the object. The 
routine getshapes then takes Vg as the size of the gaussian 
weight function used to measure quadrupole moments lij 
of the object about its center of light. The ellipticity of the 
object is then defined as d = {hi — I22, 2/i2}/(/ii + 122)- 

The next step in the KSB algorithm is to correct for 
the anisotropy of the PSF. The corrected ellipticity of a 
galaxy e^^j-j.^^^^^^ is related to the observed smeared ellipticity 
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where the eUipticities are understood to denote the rele- 
vant 2-component spinor e^, and p is a measure of PSF 
anisotropy. The tensor Pg^n is the smear polarisabihty, a 
2x2 matrix with components involving higher moments of 
surface brightness. Since for stars e*^^^^^^^^ = Q, p can be 
measured using 

P ^^sm) ^smeared- (3) 

The lensing shear takes effect before the circular smear- 
ing of the PSF. Luppino and Kaiser (1997) showed that the 
pre-smear shear 7 averaged over a field can be recovered 
using 

T ^1 ^corrected (^) 

where 

p* 

P — P9 _ sfe p9 irC\ 

Here, P^i_^ is the shear polarisabihty tensor for the galaxy 
involving other higher order moments of the galaxy image. 
The quantities P*^ and P*„ are the shear and smear polar- 
isabilities calculated for a star interpolated to the position 
of the galaxy in question. With the smear and shear po- 
larisabilities calculated by imcat, we can therefore find an 
estimator for the mean shear in a given cell. 



elements are small and the diagonal elements are equal 
within the measurement noise (typical Ps*m, 11,22 = 0.10, 

^s*m, 12,21 < 5 X 10 , P*hfiO,l\ ~ 1-1) Psh,\2,2l < 0.01). 

With this simplification, we calculate P^, according to 
equation Pry is typically a noisy quantity, so we fit it as 
a function of rg. We choose to treat P-y as a scalar, since the 
information it carries is primarily a correction for the size of 
a given galaxy, regardless of its ellipticity or orientation. We 
thus plot P^^ and P^^ together against rg, and fit a cubic 
to the combined points. Moreover, since P^ is unreliable for 
objects with Vg measured to be less that rg, we remove all 
such objects from our prospective galaxy catalogue. Finally, 
we calculate a shear measure for each galaxy as in Equa- 
tion (^, where the P^ is the fitted value for the galaxy in 
question. 

Because of pixel noise, a few galaxies yield extreme, 
unphysical, shears 7. To prevent these from unnecessarily 
dominating the analysis, we have removed galaxies with 7 > 
2. 

This entire procedure provides us with an estimator 
of the shear 7 for each galaxy. We can also calculate the 
mean shear 7 ~ il) in ^ cell and its associated error 
a [7] ~ a-['y]/\/7^, where A'' is the number of galaxies in a 
cell. 



4.2 Specific Implementation 

Firstly, we need to remove noisy detections. We applied a 
size limit Vg > 1.0 to reject extraneous detections of very 
small objects claimed by im cat . We also applied a signal- 
to-noise v > 15.0 limit (see ^5.2| for justification of this ap- 
parently very conservative cut). To reduce the noise in our 
measurement, we also remove highly elliptical objects with 
e > 0.5. 

Stars were identified using the non-saturated stellar lo- 
cus on the magnitude-rji plane (see figure 11 in BRE), 
typically with R ~ 19-22. In the data, the stellar elliptic- 
ity is a smooth function of position on the field. We thus 
adopted an iterative interpolation scheme to model this vari- 
ation. Specifically, we first fitted a 2-D cubic to the mea- 
sured stellar eUipticities, plotted the residual eUipticities 
gres _ g* _ g/it re-fitted after the removal of extreme 
outliers (caused by galaxy contamination, blended images 
and noise). The stellar ellipticity was kept constant in the 
simulations, but we nevertheless fit the 2-D cubic for cor- 
rection, as a means of retaining potential systematic efi'ects 
induced in this step. 

In order to correct galaxies for anisotropic smear, we not 
only need the fitted stellar ellipticity field, but also the four 
component stellar smear and shear polarisabilities as a func- 
tion of position. Here a 2-D cubic is fit for each component of 
Psm and P*f^ . Galaxies are then chosen from the magnitude- 
r;j diagram by removing the stellar locus and objects with 
V < 15, rg < 1, e > 0.5, as described above. From our fitted 
stellar models, we then calculate e*, P^^ and P*^ at each 
galaxy position, and correct the galaxies for the anisotropic 
PSF using equation (^). As a result, we obtain 6^^,.,.^^^^^ for 
all selected galaxies in each cell. 

We then calculate P-y for the galaxies. We opt to treat 
P*^ and P*„ as scalars equal to half the trace of the re- 
spective matrices. This is allowable, since the non-diagonal 



5 RESULTS 

5.1 Ellipticity Distribution 

We first compare the distribution of eUipticities and shear 
estimators within a field. As we noted above, the uncorrected 
ellipticity distribution of the simulated objects is very sim- 
ilar to that of the data (Figure ^). The distribution of the 
shear estimators 71 after all corrections is shown on Fig- 
ure for each case. The agreement is good, showing that 
the simulations faithfully reproduce the shape statistics of 
the data, the central concern for weak lensing. 

In general, the variance of the shear estimators results 
from several effects: the intrinsic ellipticity dispersion, pix- 
elisation, and pixel noise. The latter effect is enhanced by 
the correction for the isotropic smearing. The ellipiticity and 
shear rms dispersion, ae = {e^)2 and a-y = {7^)2, at differ- 
ent stages of the correction algorithm are listed in Table |l| 
The ellipticity dispersion observed in the simulated image 
is reduced by the PSF smearing as compared to the input 
ellipticity dispersion. (Pixel noise and pixelisation tends to 
increase the observed dispersion, but the smearing domi- 
nates). The smearing is corrected for by the KSB method, 
leading to a re-increased dispersion in the corrected shear 
estimator. 

We can obtain an estimate of the relative contribution 
of pixel noise and intrinsic dispersion using these results. 
In the absence of weighting and smearing, the ellipticity is 
related to the shear by e = (77, where g = (2 — (e^)) (see 
Eq. [0] and Rhodes et al. 2000). For the simulations, the 
input ellipticity dispersion is (Je — 0.47, yielding g ~ 1.8. As 
a result, the shear rms produced by the intrinsic dispersion 
alone is ^ cr^^P^Vg ~ 0.26. The fact that this value 

is close to the total shear disperion a^y ~ 0.31 observed in the 
simulations shows that the intrinsic dispersion is larger but 
comparable to that produced by pixel noise and pixelisation. 
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Table 1 . EUipticity and Shear rms dispersion at different stages pf the correction method 



dispersion 


simulation 


data 


(Te (input) 


0.47 




(Te (measured) 


0.20 


0.20 


cr-y (final) 


0.31 


0.39 



This considerations should be kept in mind in planning the 
exposure time of weak lensing surveys. We will study the 
i mpa ct of worsened seeing and larger pixel sizes in § 5.4 and 



5.6 



5.2 Test of the anisotropic correction 

Before we discuss the reclamation of shear by the method, 
we address the existence of a remaining systematic effect. In 
BRE, we found that a signal/noise cut oi v > b (as opposed 
to our conservative u > 15) reveals a strong anti-correlation 
between the mean shear 7^ and the mean stellar ellipticity 
e*. Here we show that the same effect is found in the simu- 
lated data. Figure ^ shows 7^ vs. e* for 20 simulated fields, 
which exhibit similar behaviour to that found in BRE sec- 
tion 6 for real fields. To assess the significance of this effect, 
we use the correlation coefficient 

a = ^^'^'^-f^^-^'^ (6) 
cr(e*)a(7,) 

For a > 5 cut we find Ci = -0.69, C2 = -0.81 for 20 cells, 
which corresponds to a 3> 3(t effect as for the real data. A 
cut at > 15 reduces the anti-correlation to Ci = —0.32, 
C2 = —0.48. This corresponds to a 1.5-2(t effect, which is no 
longer a significant contribution to the lensing amplitude. 
This anti-correlation is thus due to an over-correction of the 
PSF for small galaxies (in Eq. |^). This is likely to arise from 
the fact that, because of noise, the observed radius (and thus 
(^sm)~^) of faint galaxies is slightly smaller than that of the 
bright stars used to measure the PSF. Note that figure ^ 
is very similar to the equivalent figure in BRE section 6 for 
real data. This again confirms the validity of the simulations 
and their use in testing systematic effects, and verifies that 
the low-level anti-correlation found in real data stems from 
a reproducible problem with the current correction method. 



5.3 Shear Recovery 

We now wish to observe how the output shear derived by 
the KSB algorithm compares to the shear input. In order to 
test this, we ran a set of simulations with 5% rms shear for 
20 fields. More precisely, for each field we drew a uniform 
shear from a Gaussian probability distribution, with stan- 
dard deviation equal to a shear of 0.05. A similar set of 30 
fields were simulated with 1.5% rms shear. 

Our results for the 5% simulations are shown in fig- 
ure ^ (see also section 7 in BRE for a summary). The figure 
shows that the output shear is clearly linearly related to 
the input shear, with a slope close to 1. As a quantitative 
test, we apply a linear regression fit to both components of 
the shear, combined and separately. For the combined com- 
ponents we obtain 7°"' = 0.0007 + 0.847^", with standard 
errors on the coefficients of 0.001 and 0.04, respectively. For 




0.2 




0.2 



Figure 6. The anti-correlation of e* and 7^ plotted for 20 simu- 
lated fields, where i=l,2 have been superposed, for (top) a. v > b 
cut, and (bottom) a, v > lb cut. Note the trend for v > b. 



the individual components we obtain 7°"' = 0.002 + 0.907^° 
with errors (.001, .05) and 7!"* = 0.0001 + 0.767^" with er- 
rors (.001, .04). For the 1.5% simulations we similarly obtain 
consistent results, namely 7°"* = 0.0001 + 0.797^° for com- 
bined components with respective standard errors of 0.001 
and 0.091. 

We see that the imcat measure of shear is symmetrical 
about zero, but is measuring a slightly smaller shear signal 
than the input shear. In similar conditions, we should there- 
fore adjust our shear measures by dividing 71 by 0.9 ± 0.05 
and 72 by 0.76 ± 0.04 when using this KSB implementa- 
tion. A full discussion of the recovery of rms shears using 
an extensive statistical analysis can be found in BRE, in- 
cluding discussion of the recovery of rms shears from sets of 
simulated fields. 



5.4 Effect of Seeing 

Of great practical interest is the dependence of the sensitiv- 
ity of weak lensing measurements on seeing. To study this 
dependence, we ran several simulations with the same object 
catalogue, but with different seeing values, for a set expo- 
sure time of 1 hour. For each simulated 8' x 8' simulated 
field, we computed the rms noise cTnoiso = o-y/^/N , where 
is the rms of shear measures in a single field, and A'" is the 
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Figure 7. compared with 7°"' for simulated data sheared by 
5% rms shear; top panel shows 71 input and output, while bottom 
panel shows 72 component. The dashed line shows the 7^" = 7°"' 
relation; the solid lines shows the best fits, 7J" = 0.907°"' 
7^" = 0.767°"*. 



number of usable galaxies in the field. The quantity (Tnoisc 
is a measure of the uncertainty for measuring the average 
shear in the field. 

The results are shown in table ^ and figure |^. As can be 
seen in the figure, the seeing degrades the uncertainty almost 
linearly. Interestingly, the loss of sensitivity comes primar- 
ily from the loss in the number A'^ of usable galaxies, with 
no strong increase observed in (see table El) . One might 
suppose that this degradation could be countered by longer 
integrations on the field, to regain the number counts diluted 
by the larger isotropic smear. However, besides the integra- 
tion time increase being conside rabl e for such a reclamation 
of number density (see section |5.5| ), many of the regained 
galaxies will still need to be excluded as their shape infor- 
mation has been erased by a kernel significantly larger than 
their intrinsic radius. The noise could perhaps be reduced 
by improved shear-measurement methods, which reduce the 
cuts we have to make on small galaxies. 

Note that, for worse seeing cases, the usable galaxies 
will be on average brighter and larger and will thus have a 
lower median redshift. This will tend to degrade the lensing 
signal further. For a cluster normalised CDM model, the 
shear rms from lensing in an 8'x8' cell is aions — 0.012z'^^ 
(BRE). The median redshift Zm of the galaxies is derived 
from the median 7?-magnitude using the results of Cohen et 
al (2000). The resulting lensing rms a-i^ns is also plotted as 
a function of seeing in figure ^, and the signal-to-noise ratio 
for a single (8'x8') cell, S/N — criens/o"noisc is listed in table 
^ We find that the reduction of aiens with poorer seeing is 
rather weak. Thus the reduction of signal-to-noise for shear 
measurement is again dominated by the decrease in A'^. 
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Figure 8. (Tnoisc (solid) and cicns (dashed) as a function of seeing 
FWHM for a set of simulations with 1 hour integration time. Note 
the steady degradation of sensitivity with seeing. 



5.5 Effect of Integration Time 

To optimise weak lensing surveys, one needs to compro- 
mise between depth and width. To help in this optimisation, 
we produced several simulated images for different exposure 
times, while keeping the seeing at 0.8". Table ^ shows the 
quantities discussed in the previous section for different ex- 
posure times relevant for ground-based observations. The 
noise and lensing rms are plotted on figure ^. The depen- 
dence of these quantities on exposure time is rather weak. 
This is due to the fact that the fainter galaxies which can be 
detected with deeper exposures are too small to be resolved 
in the presence of seeing, and must therefore be mostly dis- 
carded. Moreover, since intrinsic ellipticities dominate the 
total ellipticity dispersion in this regime (see §5.1), the re- 



duced pixel noise of deeper images does not substantially re- 
duce (Tnoisc- As a result, the signal-to-noise S/N to measure 
lensing is only moderately improved for longer exposures. 



5.6 Effect of Pixelisation 

Another test of practical interest is the dependence of the 
sensitivity to lensing on pixel size. This is important for the 
design of future dedicated instruments (eg. Kaiser et al 2000, 
Tyson et al 2000). To study this dependence in the context 
of current, ground-based observations, we again produced 
several simulated fields from the same object catalog, keep- 
ing a seeing FWHM of 0.6" and exposure time of 1 hour, 
but with different pixel sizes, ranging from 0.1" to 1.0". 

The results are listed in Table ^ and Figure We find 
that the noise is quite stable for pixel scales smaller than 
the seeing radius. Increased oversampling of the PSF does 
not improve the noise properties to any great degree. How- 
ever, as the pixel scale increases above the seeing FWHM, 
the method fails quickly. The stellar locus on a magnitude- 
radius plot then approaches the galaxy locus even at bright 
magnitudes, making star selection for anisotropic smear cor- 
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Table 2. Shear sensitivity as a function of seeing, in {8'x8') cells. 



Seeing (") 


Ug (arcmin 




'^noisc 


Median R 


Median z 


'-''lens 


S/N 


0.4 


29.8 


0.43 


0.0097 


24.1 


0.8 


0.0096 


1 n 


0.8 


18.0 


0.44 


0.0130 


23.4 


0.8 


0.0096 


0.7 


1.2 


11.8 


0.46 


0.0168 


22.8 


0.7 


0.0086 


0.5 


1.6 


7.5 


0.49 


0.0225 


22.2 


0.6 


0.0076 


0.3 


2.0 


5.2 


0.50 


0.0275 


21.9 


0.6 


0.0076 




Table 3. Shear sensitivity as 


a function of integration time, for (8'x8') cells. 




time (s) 


Ug (arcmin ^) 




^ noise 


Median R 


Median z 


^lens 


S/N 


1800 


11.3 
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0.0160 
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0.0086 
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0.42 


0.0117 


23.6 


0.8 


0.0096 


0.8 
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22.0 
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0.0111 
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Figure 9. Effect of integration time on the noise (solid) and 
lensing (dashed) rms. 



Figure 10. Effect of pixel size on the noise (solid) and lensing 
(dashed) rms. 



rection very difficult. Moreover, a pixel scale of, say, 0.8" 
with a seeing FWHM of 0.6" are both conspiring together 
to remove shape information grievously from galaxies be- 
yond R ~ 22. As a result, (Tnoiso rapidly grows. Extreme 
oversampling of the PSF appears to be inefficient, while 
under-sampling is very detrimental for typical ground-based 
seeing. Under-sampling is less of a problem for space-based 
data, however, since the typical pixel scale and PSF FWHM 
(say 0.1") are much smaller than the typical galaxy radius 
at magnitudes of interest (eg. Rhodes, Refregier & Groth 
2000). 



5.7 Effect of Overlapping Isophotes 

Spurious lensing signals could also be produced on small 
scales by overlapping isophotes of neighboring galaxies. Van 
Waerbeke et al. (2000) suggested that this effect could ex- 



plain the excess small scale power observed in their cosmic 
shear surveys. To understand this effect, consider two galax- 
ies which, for simplicity, are assumed to be circular. If the 
galaxies are separated by an angular distance comparable 
to their angular sizes, their combined isophotes will have a 
dumb-bell shape rather than be the simple sum of two dis- 
joint circular isophotes. As a result, their ellipticity will tend 
to be aligned along the separation axis, leading to a spurious 
ellipticity correlation. 

To test the impact of this effect, we performed 10 sim- 
ulations without lensing shear (with seeing 0.8" and expo- 
sure time of 1 hour). For these simulations, we measured the 
shear correlation functions defined as (eg. Kamionkowski et 
al. 1998) 

c^Ae) = {1^{QhA0)), (7) 

where i and j run from 1 to 2. Here, 7i(0) and 7j(^) are 
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Table 4. Shear sensitivity as a function of pixel size, for (8'x8') cells 



Pixel (") 


Ug (arcmin ^) 




•^noisc 


Median R 


Median z 


^Icns 


S/N 


0.1 


23.9 


0.43 


0.0126 


23.5 


0.8 


0.0096 


0.8 


0.2 


23.4 


0.33 


0.0099 


23.5 


0.8 


0.0096 


1.0 


0.4 


21.7 


0.35 


0.0107 


23.4 


0.8 


0.0096 


0.9 


0.6 


19.5 


0.35 


0.0114 


23.4 


0.8 


0.0096 


0.8 


0.8 


12.1 


1.02 


0.0367 


23.1 


0.7 


0.0086 


0.2 


1.0 


13.7 


0.66 


0.0252 


23.0 


0.7 


0.0086 


0.3 



the shear estimates of each member of a galaxy pair with a 
separation 0. These ellipticities are measured in a coordinate 
system whose x-axis is along the separation angle of the 
galaxy pair (see Heavens, Refregier & Heymans (2000) for 
an illustration) . After measuring these correlation functions 
for each of the 10 simulations separately, we computed the 
mean and error in the mean over all simulations. 

The resulting correlation functions are shown on Figure 
0. As a comparison, the correlation function expected from 
lensing for a cluster-normalised KCDM model is also shown 
(see ibid for details of the calculation). In the null simula- 
tion, we weakly detect a correlation on scales smaller than 
about 1'. The amplitude of the overlapping isophote effect 
is Ci ~ 10~* corresponding to an rms shear of about 1%. 
As can be seen on the figure, this is smaller, but compara- 
ble to the lensing signal expected on these scales. The exact 
amplitude of the overlapping isophote effect will depend on 
the precise conditions of the observation (or simulation). It 
is nevertheless likely that the excess power observed by van 
Waerbeke et al. (2000) on small scales [6 ^ 10") is indeed 
due to this effect. Apart from residual systematic effects, an- 
other explanation for this excess power could be the intrinsic 
alignment of galaxies. Theoretical studies however indicate 
that this effect is small for a survey of this depth (Heavens, 
Refregier & Heymans 2000; Metzler & Croft 2000; see also 
Catelan, Kamionkowski & Blandford 2000, and Pen, Lee & 
Seljak 2000). 



6 CONCLUSIONS 

We have tested our shear measurement method using numer- 
ical simulations of artificial images. The object catalogues 
were created by generating realisations of the HST Groth 
Strip; the resulting artificial images include observational ef- 
fects such as noise, seeing, and anisotropic PSF. We compare 
our realised catalogues to those observed with WHT and find 
good statistical agreement. We used these simulations to test 
the KSB shear measurement method. Overall, we find that 
this method is rather accurate, but with several provisos: we 
find a residual anti-correlation between the PSF ellipticity 
and the corrected ellipticities of faint galaxies. This effect 
can be made negligible if faint galaxies (with S/N ^ 15) are 
removed from the catalogue. We also find that the recov- 
ered shear is linearly related to the input shear, but with 
a coefficient of about 0.8 which must be used to calibrate 
the final shear. With these precautions, the KSB method is 
sufficient for the current weak lensing surveys. However, the 
method is neither optimal nor necessarily extendable to su- 
perior observing conditions, ft should therefore be replaced 
with more accurate methods such as that of Kaiser (1999b), 
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Figure 11. Mean correlation functions Ci and C2 for 10 null 
simulation fields. Points show the mean of 10 correlation func- 
tions, with associated uncertainty. The dash-dotted curves are the 
correlation functions for the expected lensing signal for a cluster- 
normalised ACDM cosmology. 



Rhodes, Refregier & Groth (2000), and Kuijken (1999), in 
future, more sensitive surveys. 

We also used our simulations to study the effect of see- 
ing, exposure time and pixelisation on the sensitivity to the 
shear. We found that increased seeing FWHM increases the 
noise almost linearly, with the primary loss being the de- 
creased number of usable galaxies. In the seeing-dominated 
regime, the sensitivity to shear is only weakly dependent 
on exposure time. As long as this regime holds, it is there- 
fore more efficient to tend towards wide rather than deep 
weak lensing surveys. Increased pixel scale hardly affects the 
sensitivity until the pixel scale is comparable to the seeing 
FWHM, at which point the method fails for typical ground- 
based seeing. Thus, extreme oversampling of the PSF does 
not seem to be necessary. 

We also tested the claim by van Waerbeke et al. that 
spurious shear signals on small scales (6^ ^ 10") could be 
produced by overlapping isophotes of neighboring galaxies. 
Using simulated images without input shear, we weakly de- 
tect this effect on scales 9^1'. The rms amplitude of this 
effect is of about 1%, which is smaller but comparable to 
that expected for lensing. Overlapping isophotes are thus 
likely to explain the excess power found on small scale by 
this group. 
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